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We present a calculation of open charm quark production in diffractive deep inelastic electron- 
proton scattering in a perturbative QCD framework. The cross section is proportional to the 
square of the gluon density and explicitly displays breaking of Regge factorization. Jet cross 
sections as well as the charm contribution to the diffractive structure function are calculated. As 
a consequence of the steep rise of the gluon density at small x the charm contribution to F2 
rises with decreasing xp. 
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1. ) In the context of the discussion of rapidity gap events in deep inelastic electron-proton scattering observed 
at HERA a subclass of diffractive events in which a large mass scale appears in the diffractively produced 
hadronic final state has created particular interest. Representatives of this type of events are the diffractive 
production of heavy vector mesons diffractive jet production ||^, |4| ||, ^ and diffractive production of 
open charm ||] . Due to the presence of the large mass scale these processes offer the possibility to apply 
and test perturbative QCD in the setting of diffractive scattering. The common feature of the above cited 
processes is the dependence of the cross section on the square of the gluon density of the proton. Because of 
this strong sensitivity these events have been considered as a possible probe of the gluon density. 

In the above list the process of diffractive open charm production is particularly promising. Compared to 
heavy vector meson production it does not depend on a meson wave function which is poorly determined 
from the theoretical side. Compared to jet production charm production does not require a large transverse 
momentum which in turn leads to a strongly suppressed event rate. 

In this letter we generalize our preceding calculations on jet production ^ in DIS diffractive dissociation 
to finite quark mass. This allows, in particular, the investigation of open charm production. Our calculation 
is based on an analytical expression for the unintegrated gluon structure function which enables us to take 
into account a subset of subleading corrections which turn out to be numerically important Q . As the new 
contribution of the present work we calculate jet cross sections for charm quarks and compare with massless 
flavours. In addition we compare the magnitude of the jet cross section calculated in our model with the 
predictions of a model which is based on nonperturbative two-gluon exchange [Q . Furthermore our formulae 
for the elecron-proton cross section include the dependence of the cross section on the angle of the jet plane 
relative to the electron plane. For the case of the large charm mass we can extend our expressions to low 
transverse momenta and can even integrate the transverse momentum to obtain the charm contribution of 
the diffractive structure function. We discuss the xp , /3 and dependence of the charm contribution to 
F^^-^\ In this part of our analysis we obtain results similar to the ones presented in Q and 

2. ) The kinematics of the process is well-known and we only give the key ingredients. To calculate first the 
hadronic tensor H^^ consider the photon-proton subprocess 7* (q) + P{p) — > c{k) + c{q — fc -I- A) -I- P{p — A). 
We define the momenta of the particles as indicated in figure |l|. The calculation is performed for zero 




Figure 1: Diffractive production of a quark- antiquark pair in photon- proton scattering. 



momentum transfer t = A^ = 0. For the momenta k and A we use a Sudakov decomposition w. r. t. the 
light cone momenta p and q' = q + xp {x = /{2pq)) 

k = aq' + l3p + k, (1) 
A = aAq' + XPP + A. (2) 

Using the mass shell conditions for the outgoing particles one can show that we have aA = 0(A^/(2pq')). 
Assuming <^ pq' and hence neglecting aA and A we can cast the phase space in the form 

dT = — ^ rlM"^ dt dk^ (3) 

8pq L . k2+m2 
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with the charm quark mass and being the invariant mass of the cc pair which is related to the hght 
cone momentum fraction a through the relation 

- (4) 
a(l — a) 

Energy- momentum conservation then leads to the phase space restriction > 4(k^ + ™c)- Furthermore 
the longitudinal momentum fraction xp transferred from the proton to the cc pair is fixed as 

W'^ = {p + qf' is the cms-energy of the photon-proton system. Another often used variable is (3 defined as 

from which follows (3 = x/xp. 

Now we have to specify the coupling of the cc pair to the proton. As the simplest model for an interaction 
in which no color is transferred from the proton to the cc pair we take two-gluon exchange. Since the 
charm quark is sufficiently heavy we treat both gluons perturbatively. We then make use of high-energy 
factorization |^ to express the amplitude of the photon-proton subprocess in terms of the unintegrated gluon 
density of the proton (fig. |^) 

H^, = (^j dHC^{\-]^,Q\M^) Tg{xp,\^)^ (^j dH'Cu{\'-\,,Q\M^)TG{xp,\''')^ . (7) 

This factorization is valid in the leading-log(l/a;jp) approximation in which the imaginary part of the di- 
agrams in fig. 1^ contributes. In this approximation the difference of the longitudinal momenta of the two 
gluons is neglected. It is therefore legitimate to use the same diagonal gluon density !Fq{xp,\^) which ap- 
pears in inclusive DIS. The coefficient function represents the upper part of the diagrams in fig. S with 




Figure 2: Representation of the amplitude in terms of the unintegrated gluon structure function. The dots 
represent three other diagrams which are generated by attaching the gluons to the quark lines in all possible 
ways. 



the gluons coupled to the quarks in all possible ways. The calculation of the respective diagrams is straight- 
forward if one uses the fact that the gluon polarization proportional to p'^ gives the dominant contribution 
in the small-x limit |lO[ [ill 0. 

In order to obtain the electron-proton cross section we have to contract Hf^i, with the usual lepton tensor. 
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This yields the foUowing expression 



da. 



2yQ^Tr^ 
+ (1-2/) 
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\t=o 
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|t=o 



(8) 



The indices T and L refer to the contributions of transversely and longitudinally polarized photons. The term 
with index / is the interference contribution. The angle (j) is defined in the two-dimensional transverse plane 
through the scalar product e-ki = |e| |k| cos 4>. Here e is the transverse component of the electron momentum 
= l/yq'^ + x{l — y)/yp^ + and ki is the transverse momentum of the charm (anti)quark which points 
into the proton hemisphere (in the cc cm-system) . Note the negative sign of the cos 2(/)-term which is a 
characteristic feature of the two-gluon exchange model. In photon-gluon fusion (one-gluon exchange) this 
term comes with positive sign 

The expressions for the subprocess cross sections read 
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The factor Cc denotes the charge of the charm quark. The essential dynamics is contained in the universal 
functions II , It for which we have the expressions 
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with Pk.i being defined as 

k2 



Pki = (1' 



2\\2 



M^)) 



4k^(^0^ 



+ ^(0' + M')) 



(14) 



(15) 



These expressions simplify considerably in the case of massless flavors (m^ =0) || . To proceed one has 
to use an explicit model for the unintegrated gluon structure function Tq- Within the leading-log(l/xip) 
approximation it is consistent to express as the solution of the BFKL equation |Q . The 1- integration in 
eqs. (^ 14) can then be performed analytically. The main outcome of this calculation is the determination 
of the relevant scale A2 of the process which is found to be 



A2 = (k2+m2) 



(16) 
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If we assume that this effective scale is large (compared to the inverse size of the proton) we can evaluate 
the integrals in eqs. (pUlS) in the leading-log(A^) approximation. Using the relation 



dl^TaixpA^) = xpGixp.A'') (17) 



we can express the functions Il,It within this approximation in terms of the gluon density X]pG{x]p, A^) 
The result reads 

2^n2 f_1 1,2 , ^2 



= k^(rKF"^^("^'T^)- 

Here we have introduced, for brevity, the scaling variable £, — {(3 + TO^/k^)/(l — (3). 

Inserting these expressions into the cross section formulae we obtain the cross section for diffractive cc 
production in DIS in the double leading logarithmic approximation (DLA). In a certain class of corrections 
to the DLA was obtained in the massless case which turned out to be numerically important in the small-/? 
region. The same type of correction can be calculated in the massive case and the result for the correction 
terms reads 

.(e) _ {^'+m))Q' 1 , ? \ d , 



= k^AP (l+a^ (^2 + (l-01og^j^XfG(x^,k-(l+0) (20) 
= 2(^^(l-e-2eiog^)^xrG(xf,k2(l + C)) (21) 



These terms have to be added to the respective contributions ( p^ ) and (p^. In the next section we turn to 
the numerical evaluation of our formulae. 

3.) As to the numerical evaluation we start with some preliminary remarks. In all subsequent calculations 
we have integrated over the momentum transfer t after multiplication of the results of 2.) with the elastic 
proton form factor 

^P(---^) = "r'-'r" fl - TT^) C"^- (22) 



GcV^ 



0.7 GeV' 



For the charm quark mass we have taken the value mc = 1.5 GeV and we have used the one- loop formula for 
as evaluated at the scale (k^ -I- m^)/(l — /?). In a preceding publication we have discussed whether the 
leading order or the next-to-leading order gluon density should be used for a numerical prediction. We have 
presented some evidence that the NLO density is appropriate and consequently we perform all calculations 



with the GRV |T5| NLO gluon density. 

First we present results for jet cross sections, i. e. we calculate the total electron proton cross section with 
a lower cutoff kQ > 2 GeV^ imposed on the transverse momentum integration. The other variables are 
integrated in a region which is adapted to the HERA kinematics, namely we choose the cuts > 10 GeV^, 
xp < 10"^ and 50 GeV < W < 220 GeV. The results for the integrated cross section are listed in table |^. 
For comparison we quote in table ^ also the results which are obtained for three massless quarks . 
It can be clearly seen that the relative contribution of the charm quarks compared to the massless flavors 
is dependent on the lower momentum cutoff kg. The ratio r = Cce/ X]/=m d s "^/Z rises from r — 0.25 for 
kg = 2GeV^ to r = 0.53 for kg = 8GcV^. This is easily understood since in the limit m^/kg the ratio 
r should approach 2/3 which follows from counting the electromagnetic charges. It is also noticeable that 
the decrease of the transverse cross section with increasing kg is weaker in the case of charm production 
compared to the massless case. Stated differently, the k^-spectrum of the transverse contribution is flatter in 
the massive case. As a consequence the ratio of the longitudinal to the transverse contribution is decreasing 
with increasing kg for charm production whereas for the massless flavors it is approximately constant. 
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= 2 GcV^ 


k^ = 4GcV^ 


ki^ = 8 GeV^ 


DLA + corrections 


< 


27 


11 


3.4 




2.3 


0.4 


0.1 




29.3 


11.4 


3.5 



Table 1: Results for total eP-cross sections (in pbarn) of diffractive dijet production for charm quarks. 





k^ = 2 GcV^ 


kj^ = 4GeV^ 


kj^ = 8 GcV^ 


DLA + corrections 




108 


30 


6 




9 


2 


0.6 




117 


32 


6.6 



Table 2: Results for total eP-cross sections (in pbarn) of diffractive dijet production for three massless 
flavors. 



The absolute numbers have been compared to the ones computed in a model based on nonperturbative two- 
gluon exchange (soft pomeron model) already in |jl^. The soft pomeron model predicts cross sections 
which are reduced by a factor of approximately 1/2 compared to our calculation. This difference is mainly 
due to the rise of the cross section at small xp m. our calculation according to the rise of the function 
\xipG{x]p, . The soft pomeron model does not contain such a pronounced rise at small xp. 
So far we have discussed jet cross sections. Since the large charm quark mass TOc provides a hard scale 
independent of the value of the transverse momentum k^ we can also extend the perturbative calculation 
to low values of k^ . We can even integrate over the transverse momentum which gives us the possibility to 



calculate the charm contribution to the diffractive structure function P 



D(3) 



in our model. The latter has 



been introduced to describe diffractive events (events with a rapidity gap) in DIS. In terms of Pj 
diffractive electron proton cross section reads 



the 
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where the longitudinal contribution has been neglected. With this definition of daf^ 
from the photon-proton cross section in the following way 
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can be obtained 



(24) 



where the cross sections in brackets are obtained from multiplying the expression for t = {) in eqs. (|^) and 
( p^ ) with the proton form factor Fp{xp,t) and the factor 27r (integration of 0). In fig. |^ we display P^*-^-* 
as a function oi xp. First of all, these curves demonstrate the rapid increase of the diffractive structure 
function at low values of xp . The fact that the shapes of the curves are not identical for different [3 indicates 
the (weak) breaking of Regge factorization. Variation of (3 leads to a variation of the xjp-dependence. More 
precisely, a lowering of /3 leads to an increase of the slope of the curve in accordance with the fact that the 
relevant scale of the gluon density is proportional to 1/(1 — P). For smaller (left hand side) this effect 
is very small. By comparing the results for — 20 GeV^ and — 50 GeV^ one observes that a variation 

This is demonstrated 



D{3) 



of does not lead to a significant change of the absolute normalization of Fr 
explicitly in fig. where the Q'^-dependence of F.^^^\ respectively its transverse and longitudinal part, is 
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Figure 3: The Xip-dependence of the difFractive structure function for ~ 20GeV^, — 50GeV^ and 
/3 = 1/3, (soUd hne) and P = 2/3 (dashed hne). 
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Figure 4: The Q^-dependence of the diffractive structure function for (3 — 2/3 and xp = 10^^. The dotted 
hne is the transverse contribution the dashed hne is the longitudinal part and the solid lines represent the 
sum. 



displayed. One observes a weak logarithmic increase of the transverse contribution whereas the longitudinal 
part decreases rapidly with Q^. The longitudinal part, in other words, is a higher twist contribution. 
In fig. H we show /3-spectra for two different values of Q^. The most noticeable feature here is of course 
the charm threshold at the right end of each curve. One sees that near this threshold (for large /3) the 
longitudinal contribution becomes comparable to or even larger than the transverse contribution. In the 
limit /3 — > both the transverse and the longitudinal part tend to zero. 

We can compare our results with experimental data on inclusive diffractive DIS to roughly estimate the rel- 
ative contribution of the cc-production. The ZEUS collaboration gives a value of 30-40 for F^^^\l3, Q^,X]p) 
at xp = IQ-^, Q2 ^ 16GeV^ and /? = 0.65 0. From the right hand side of fig. | we find a value of 



6 




Figure 5: The /3-dependence of the diffractive structure function for — 20 GeV^ and ~ 50 GeV^ and 
xp = 10~^. The dotted Hne is the transverse contribution the dashed hne is the longitudinal part and the 
solid line represents the sum. 



^2^^^^ ~ 3.25 for xip = 10-3, ^ 20GeV^ and /3 ~ 0.65. We conclude that approximately 10% of the 
diffractive events observed at HERA in the corresponding kinematical region are due to cc-production. A 
similar estimate has been given in [Q. One should emphasize that the relative contribution of cc is a func- 
tion of xp due to the strong rise of the cc-cross section at low xp. Since the cross section of the diffractive 
production of massless flavors is not expected to increase equally fast at low xp due to the dominance of 
low momentum scales, the cc contribution is enhanced at small xp. The relative contribution of 10% seems 
to be large enough to expect the charm threshold to be observable in the experimenta data. 
4.) To summarize, we have calculated the cross section for the diffractive production of an open cc-pair in 
DIS. We have assumed two-gluon exchange and have treated both gluons perturbatively due to the large 
charm quark mass. In the double logarithmic approximation the cross section is proportional to the squared 
gluon density \xpG{xp, A'^)\'^ of the proton with the scale = (k^ -I- m^)/(l — We have first calculated 
jet cross sections where is kept large and found that the charm contribution rises relative to massless 
flavours with increasing k^. For very high k^ the relative charm contribution follows from the electromag- 
netic charge factors. Since is sufficiently large the transverse momentum k^ can be integrated and the 
charm contribution to the diffractive structure function can be calculated in our model. As a function of Xip 
the diffractive structure function rises steeply at small x p due to the rise of the gluon density. Since (3 enters 
explicitly the scale of the latter our formulae display the breaking of Regge factorization. As a function of 
the charm contribution to the diffractive structure function is approximately constant, i. e. it displays 
a leading-twist behavior. The most remarkable property of the /3-spectrum is the charm threshold at large 
(3. At small xp where the charm contribution is relatively largest this threshold should be observable in the 
data. At small /? both the transverse and the longitudinbal cross section tend to zero. 

As far as cc-production is concerned our calculation is complete in the given approximation. The data on 
diffractive dissociation, however, show a constant behavior in the limit f3 —>■ 0. This shows that our model 
is not adequate to describe diffractive dissociation in the complete /3-range. As an important correction the 
production of additional gluons has to be considered. Taking one additional gluon into account, one obtains 
a constant cross section at /3 = 0. Corrections of this type have been investigated in Their impact 

on diffractive charm production has been analyzed in jg]. The validity of the cited calculations has been 
limited to a restricted kinematical range. A complete calculation of the relevant corrections remains still an 
important task for future work. 



7 



Acknowledgments: I have benefited very much from discussions with Prof. J. Bartels, Dr. M. WiisthofT, Dr. 
M. Diehl and C. Ewerz. I am grateful to Prof. J. Bartels and C. Ewerz for carefully reading the manuscript. 
The support from the Studienstiftung des Deutschen Volkes is gratefully acknowledged. 



References 

[1] M.G.Ryskin, Z. Phys. C 57 (1993) 89; 

S.Brodsky, L.Frankfurt, J.F.Gunion, A.H.Mueller, M.Strikman, Phys. Rev. D 50 (1994) 3134; 
L.Frankfurt, W.Koepf, M.Strikman, Phys. Rev. D 5 4 (1996) 3194; 
M.G.Ryskin, R.G.Roberts, A.D.Martin, E.M.Levin, |hep-ph/95 11228 . 

[2] J.Bartels, H.Lotter, M.Wiisthoff, Phys. Lett. B 379 (1996) 239; 
Erratum, Phys. Lett. B 382 (1996) 449. 

J.Bartels, C.Ewerz, H.Lotter, M.Wiisthoff, Phys. Lett. B 386 (1996) 389. 
M.Diehl, Z. Phys. C 66 (1995) 181; preprint CPTH-S472-1096, |hep-ph/9610430| . 
N.N.Nikolaev, B.G.Zakharov, Phys. Lett. B 332 (1994) 177. 

M.G.Ryskin, M.Besancon, Proceedings of the HERA workshop 'Physics at HERA', Vol. 1 (edited by 
W.Buchmiiller, G.Ingelman), Hamburg 1991. 

M.Genovese, N.N.Nikolaev, B.G.Zakharov, Phys. Lett. B 378 (1996) 347. 



E.M.Levin, A.D.Martin, M.G.Ryskin, T.Teubner, hep-ph/9606443 



S.Catani, F.Hautmann, Nucl. Phys. B 427 (1994) 475; 

S.Catani, M.Ciafaloni, F.Hautmann, Phys. Lett. B 242 (1990) 97; Nucl. Phys. B 366 (1991) 135. 
A.H.Mueller, Nucl. Phys. B 335 (1990) 115. 
N.N.Nikolaev, B.G.Zakharov, Z. Phys. C 53 (1992) 331. 
E.M.Levin, M.Wiisthoff, Phys. Rev. D 50 (1994) 4306. 
L.N.Lipatov, Sou. J. Nucl. Phys. 23 (1976) 642; 

E.A.Kuraev, L.N.Lipatov, V.S.Fadin, Sov. Phys. JETP 44 (1976) 443; Sou. Phys. JETP 45 (1977) 199; 
LI.Balitskii, L.N.Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

A.Donnachie, P.V.Landshoff, Nucl. Phys. B 244 (1984) 322; Nucl. Phys. B 303 (1988) 634. 
M.Gliick, E.Reya, A.Vogt, Z. Phys. C 67 (1995) 433. 

J.Bartels, M.Diehl, C.Ewerz, H.Lotter, M.Wiisthoff, in: Future Physics at HERA, Proceedings of the 
Workshop 1995/1996 (vol. 2), edited by G.Ingelman, A.DeRoeck, R.Klanner, Hamburg 1996. 

ZEUS collaboration, Z. Phys. C 70 (1996) 391. 

J.Bartels, M.Wiisthoff, J. Phys. G 22 (1996) 929. 



8 



